An Efficient Particle Tracking Algorithm for Large-Scale Parallel Pseudo-Spectral Simulations of Turbulence

Particle tracking in large-scale numerical simulations of turbulent flows presents one of the major bottlenecks in parallel performance and scaling efficiency. Here, we describe a particle tracking algorithm for large-scale parallel pseudo-spectral simulations of turbulence which scales well up to billions of tracer particles on modern high-performance computing architectures. We summarize the standard parallel methods used to solve the fluid equations in our hybrid MPI/OpenMP implementation. As the main focus, we describe the implementation of the particle tracking algorithm and document its computational performance. To address the extensive inter-process communication required by particle tracking, we introduce a task-based approach to overlap point-to-point communications with computations, thereby enabling improved resource utilization. We characterize the computational cost as a function of the number of particles tracked and compare it with the flow field computation, showing that the cost of particle tracking is very small for typical applications

Continuation for thin film hydrodynamics and related scalar problems

This chapter illustrates how to apply continuation techniques in the analysis of a particular class of nonlinear kinetic equations that describe the time evolution through transport equations for a single scalar field like a densities or interface profiles of various types. We first systematically introduce these equations as gradient dynamics combining mass-conserving and nonmass-conserving fluxes followed by a discussion of nonvariational amendmends and a brief introduction to their analysis by numerical continuation. The approach is first applied to a number of common examples of variational equations, namely, Allen-Cahn- and Cahn-Hilliard-type equations including certain thin-film equations for partially wetting liquids on homogeneous and heterogeneous substrates as well as Swift-Hohenberg and Phase-Field-Crystal equations. Second we consider nonvariational examples as the Kuramoto-Sivashinsky equation, convective Allen-Cahn and Cahn-Hilliard equations and thin-film equations describing stationary sliding drops and a transversal front instability in a dip-coating. Through the different examples we illustrate how to employ the numerical tools provided by the packages auto07p and pde2path to determine steady, stationary and time-periodic solutions in one and two dimensions and the resulting bifurcation diagrams. The incorporation of boundary conditions and integral side conditions is also discussed as well as problem-specific implementation issues

Measurement of polarization transfer in the quasi-elastic 40Ca(e⃗,e′p⃗)^{40}{\rm Ca}(\vec{e},e' \vec{p}) process

Polarization transfer to a bound proton in polarized electron knock-out reactions, $\mathrm{A}(\vec{e},e^{\prime}\vec{p})$, is a powerful tool to look for in-medium modification of the bound proton. It requires comparison to calculations which consider the many-body effects accompanying the quasi-free process. We report here measured components $P_x^{\prime}$, $P_z^{\prime}$, and their ratio $P_x^{\prime}/P_z^{\prime}$, of polarization transfer to protons bound in $^{40}\mathrm{Ca}$, which is described well by the shell model and for which reliable calculations are available. While the calculations capture the essence of the data, our statistical precision allows us to observe deviations which cannot be explained by simple scaling, including by varying the proton electromagnetic form factor ratio $G_E/G_M$. We further explore the deviations of the ratio of the polarization transfer components from that of a free proton, $(P_x^{\prime}/P_z^{\prime})_{\rm A}/(P_x^{\prime}/P_z^{\prime})_{\rm H}$, and its dependence on the bound-proton virtuality

Brahma Is Required for Proper Expression of the Floral Repressor FLC in Arabidopsis

This is an open-access article distributed under the terms of the Creative Commons Attribution License.[Background]: BRAHMA (BRM) is a member of a family of ATPases of the SWI/SNF chromatin remodeling complexes from Arabidopsis. BRM has been previously shown to be crucial for vegetative and reproductive development. [Methodology/Principal Findings]: Here we carry out a detailed analysis of the flowering phenotype of brm mutant plants which reveals that, in addition to repressing the flowering promoting genes CONSTANS (CO), FLOWERING LOCUS T (FT) and SUPPRESSOR OF OVEREXPRESSION OF CO1 (SOC1), BRM also represses expression of the general flowering repressor FLOWERING LOCUS C (FLC). Thus, in brm mutant plants FLC expression is elevated, and FLC chromatin exhibits increased levels of histone H3 lysine 4 tri-methylation and decreased levels of H3 lysine 27 tri-methylation, indicating that BRM imposes a repressive chromatin configuration at the FLC locus. However, brm mutants display a normal vernalization response, indicating that BRM is not involved in vernalization-mediated FLC repression. Analysis of double mutants suggests that BRM is partially redundant with the autonomous pathway. Analysis of genetic interactions between BRM and the histone H2A.Z deposition machinery demonstrates that brm mutations overcome a requirement of H2A.Z for FLC activation suggesting that in the absence of BRM, a constitutively open chromatin conformation renders H2A.Z dispensable. [Conclusions/Significance]: BRM is critical for phase transition in Arabidopsis. Thus, BRM represses expression of the flowering promoting genes CO, FT and SOC1 and of the flowering repressor FLC. Our results indicate that BRM controls expression of FLC by creating a repressive chromatin configuration of the locus.This work was supported by Ministerio de Educacin y Ciencia (BFU2008-00238, CSD2006-00049), and by Junta de Andaluca (P06-CVI-01400) to J.C.R. and by the National Institutes of Health (grant no. 1R01GM079525), and the National Science Foundation (grant no. 0446440) to R.A. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.Peer reviewe

Two-particle correlations in azimuthal angle and pseudorapidity in inelastic p + p interactions at the CERN Super Proton Synchrotron

Results on two-particle ΔηΔϕ correlations in inelastic p + p interactions at 20, 31, 40, 80, and 158 GeV/c are presented. The measurements were performed using the large acceptance NA61/SHINE hadron spectrometer at the CERN Super Proton Synchrotron. The data show structures which can be attributed mainly to effects of resonance decays, momentum conservation, and quantum statistics. The results are compared with the Epos and UrQMD models.ISSN:1434-6044ISSN:1434-605